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ABSTRACT 

The [O III] A5007 line is commonly used as an indicator of the systemic 
redshift of AGNs. Also, recent studies have used the width of this emission line 
as a proxy for the stellar velocity dispersion in the host galaxy. This paper calls 
both of these assumptions into question by analyzing a sample of approximately 
400 AGN spectra from the first data release of the Sloan Digital Sky Survey. 
These spectra show that the low-ionization forbidden lines ([O II], [N II], [S II]) 
define a consistent redshift, but that the peak of the [O III] line is blueshifted in 
approximately half of the AGNs with respect to that redshift. For the sample 
studied here, the average shift is 40 km s _1 , with the largest shift being over 400 
km s -1 . The magnitude of this shift is found to be correlated with a number 
of properties, including the width of the [O III] line and the Eddington ratio 
(L/Lem), derived from the luminosity and width of FL3. 

Subject headings: galaxies: active — galaxies: Seyfert - quasars: emission lines 

1. Introduction 

It is often assumed that the narrow emission lines in active galactic nuclei (AGN) are 
indicative of some approximately stable, quiescent pool of gas that is being excited or illu- 
minated by the central source. Emission lines that originate in the narrow line region (NLR) 
are commonly used to obtain the systemic redshift (De Robertis 1985; Marziani et al. 1996; 
Richards et al. 2002) or as a proxy for the stellar velocity dispersion (Shields et al. 2003). 
Most often it is the [O III] A5007 line that is used, as it is usually the strongest and cleanest 
narrow line in the optical region of the spectrum. 

The validity of the assumptions that the A5007 line is a good indicator of systemic 
redshift or stellar velocity dispersion rests primarily on extrapolation of studies of low lumi- 
nosity objects, such as those by Heckman et al. (1981), Vrtilek (1985), Nelson and Whittle 
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(1996), and Nelson (2000). However, even these studies of the NLR characteristics in lower 
luminosity Seyfert galaxies suggest systematic departures from symmetry and stability are 
not uncommon. Heckman et al. (1981) and subsequent studies found (1) that the A5007 line 
tends to be asymmetric, with a sharper falloff to the red than to the blue, and (2) that the 
redshift of the A5007 line tends to be negative compared with the systemic velocity derived 
from the stellar absorption lines. 

Several recent studies have explored some of the characteristics of the [O III] A5007 line 
in AGNs. Zamanov et al. (2002) studied a sample of about 200 low redshift AGN, most of 
which were luminous enough to be considered QSOs. Under the assumption that the redshift 
derived from the broad H/? line is the systemic value, they found that the A5007 line gave a 
consistent redshift 90% of the time, but that, among the outliers, blueshifts were three times 
as common as redshifts. They investigated seven objects with the largest blueshifts, and 
found that they fell preferentially among the objects at the strong Fe II - narrow H/3 end of 
the Eigenvector 1 sequence (Boroson and Green 1992). In order to overcome the limitation 
of the assumption that H/3 gives the systemic redshift, they noted that for one of their blue 
outliers, I Zwicky 1, other observations of the host galaxy give a redshift consistent with the 
H/3 line and not the [O III] line. A consistent finding comes from the Veron-Cetty, Veron, & 
Goncalves (2001) study of narrow line Seyfert 1 galaxies (NLSls) in which it was found that 
approximately half of their sample had a broad, blueshifted [O III] A5007 line in addition to 
the narrow component assumed to be at the systemic velocity. 

Boroson (2003) and Nelson et al. (2004) studied the [O III] A5007 line in order to 
explore the utility of the width of this line as a surrogate for the stellar velocity dispersion. 
The Boroson (2003) study looked at a sample of 107 QSOs from the Sloan Digital Sky 
Survey (SDSS) Early Data Release (EDR), while the Nelson et al. (2004) study used new 
spectra of approximately 60 low-redshift PG QSOs. The technique used in both studies was 
to plot the width of the A5007 line against black hole mass derived using the H/? broad line 
FWHM and the luminosity at A5100 in the rest frame, Lsioo- The tight correlation between 
black hole mass and stellar velocity dispersion (Tremaine et al. 2002) in normal galaxies 
was the reference for the expected relationship. Both studies found that the [O III] widths 
show some correlation with black hole mass, but the scatter is large. Nelson et al. (2004) 
found that the more asymmetric [O III] lines have widths that are closer to those predicted 
by the Tremaine et al. (2002) fit. That study also found a correlation between [O III] line 
asymmetry and Eigenvector 1, in that objects with greater asymmetry were more likely to 
be at the strong Fe II end of the sequence, presumably correlating with higher h/lj Edd . 

The study reported here is an effort to explore in greater detail and more quantitatively 
the behavior of the [O III] A5007 line, by analyzing a larger sample of objects in which the 
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systemic velocity can be unambiguously determined. By combining a number of measure- 
ments of this line with other observed and derived properties, a better understanding of 
the geometric and kinematic structure in which the [0 III] emission arises may emerge. A 
cosmology with H = 71 km s _1 Mpc" 1 , Vt M = 0.27 and Q A = 0.73 is used throughout this 
paper. 

2. The Sample 

The sample of QSOs was drawn from Data Relase 1 (DR1) (Abazajian et al. 2003) of 
the Sloan Digital Sky Survey. The parent sample, which numbered about 3300, included all 
spectroscopically observed objects identified as QSOs having a redshift less than 0.5. From 
this list, those with spectroscopic signal-to-noise less than 15 in the R band (SN_R in the 
FITS header) were removed, leaving 804 objects. The redshift distribution and absolute 
magnitude distribution (in the I band) for the sample are shown in figure 1. The spectra 
were corrected to rest wavelengths, adopting the SDSS value of the redshift stored in the 
header of each object. 

These 804 spectra were inspected individually in the regions around the [O II] A3727 
doublet, the [N II] A6584 line, and the [S II] A6717 and A6731 lines. When visible, the 
wavelength of each of these lines was measured, by fitting a gaussian to the top half of the 
line. These wavelengths were then converted to redshifts, using the vacuum wavelengths 
3728.30A, 6585.27A, 6718.29A, and 6732.67A as the rest values for the lines. At least one 
of these lines could be measured in 654 of the 804 objects and at least two of the lines from 
different ions could be measured in 399 of the objects. 

Figure 2 shows the distributions of relative redshifts for each of the three ions. The 
histograms show the redshift of the identified line (or the average of the two [S II] lines) 
relative to the redshifts of any of the other lines that were measured. The average relative 
redshifts are 13.0, -11.2, and -3.6 km s _1 for [O II], [N II], and [S II] respectively. The 
corresponding standard deviations of the distributions are 39, 32, and 25 km s -1 . It is likely 
that the [N II] measurements are shifted a bit to the blue (negative velocities) because this 
narrow line typically lies on the red wing of the broad Ha line. Part of the apparent shift 
of the [O II] distribution to positive velocity may be due to the fact that, in most cases, its 
redshift is being calculated relative to the anomalously blue [N II] line. An additional factor 
which may also account for the greater width of the [O II] distribution is that this 'line' is 
actually a doublet with a separation between components of about 2.7A. An average value 
for the centroid of the blended components is undoubtedly not a good guess in some cases. 
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Despite the slight differences in mean and width of the distributions, they are all roughly 
symmetric, and it is likely that these lines indicate the systemic velocity of the galaxy. All 
of these lines from singly ionized species have low critical densities, and so they arise either 
in the outer parts of the narrow line region (NLR) or even further out, in the interstellar 
medium of the host galaxy itself. The 3 arcsecond diameter of the SDSS fibers corresponds 
to a physical diameter of about 10.5 kpc at the sample median redshift of 0.22. The mean 
of all visible low ionization lines in each object was adopted as its systemic redshift, and the 
spectra were adjusted to reflect rest wavelengths using these new, small offsets. 

From these corrected spectra, the wavelength of the [O III] A5007 line was measured 
in each object in a manner similar to that in which the low ionization lines were measured. 
These [O III] wavlengths were converted to a relative velocity using a vacuum rest wavelength 
of 5008. 24A for this line. Figure 3 shows the distribution of these [O III] redshifts relative to 
the systemic redshift calculated from the low ionization lines. The solid histogram represents 
all the objects for which a low ionization systemic redshift could be determined and the [O 
III] A5007 wavelength could be measured (638 objects; the total [O III] sample). The dotted 
histogram represents only the objects for which low ionization redshifts could be determined 
from two or more ions (399 objects; the high confidence [O III] sample). The bottom panel 
of this figure shows the sum of the redshift distributions for the three low ionization lines, 
in each case relative to the one or two other lines measured. 

The low ionization line distribution is well fit by a gaussian with a mean of zero (by 
construction) and a width (<r) of 35 km s -1 , a fraction of which is due to the slight shift 
between the distributions. The [O III] distribution, however, has a broad shoulder on the 
blueshifted side, with a sprinkling of objects out to -400 km s _1 . Approximately half the 
objects might be considered to be in a distribution like that for the low ionization lines with 
the other half on the blueshifted side. For the total [O III] sample, 135 (21%) of the 638 
values have positive (redshifted) relative velocities; this number is 93 (23%) of 399 values for 
the high confidence sample. Note that both the total [O III] sample and the high confidence 
[O III] sample are biased within the original sample of DR1 objects toward the low redshift 
and low luminosity end. This bias is due to the fact that the low ionization narrow lines are 
more commonly found in the lower luminosity objects. The dashed and dotted lines in the 
two panels of Figure 1 show the distributions for the total [O III] and high confidence [O III] 
samples respectively. 
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3. Properties of Objects with Different [O III] Redshifts 

In order to compare the properties of AGN with different [0 III] velocity shifts, four 
subsamples were constructed. The subsamples are composed of about a dozen objects from 
the high confidence sample with [O III] velocities from the most blueshifted extreme of the 
distribution (the VB sample), the moderately blueshifted region (the B sample), the region in 
which the [O III] velocities are in good agreement with the systemic redshift (the N sample), 
and the most redshifted extreme of the distribution (the R sample). 

For each object in these subsamples, the Fe II emission was measured and subtracted us- 
ing the template and the technique developed in Boroson and Green (1992). The continuum 
was then subtracted and the FWHM of the H/3 line and the FWHM and asymmetry index 
A20 (Whittle 1985) of the [O III] A5007 line were measured. Table 1 lists the properties of 
each object in the subsamples, including the absolute magnitude in the I band (from the I 
'PSF' magnitude listed in the SDSS DR1 database), the measured shift in the [O III] A5007 
line, the equivalent width of the Fe II multiplets between A4434 and A4684, the FWHM of 
H/9 and [O III] A5007, and the A20 value and equivalent width of the [O III] A5007 line. Also 
listed are values of log Mm and log L/L^d computed from the FWHM H/3 and the lumi- 
nosity at A5100 in the rest frame, which was calculated from the R band PSF magnitude. 
The [O III] FWHM values are corrected for the instrumental resolution, about 150 km s _1 . 
At the bottom of the table, average values of each quantity, for each of the four subsets, are 
given. 

Inititially, an attempt was made to use the tabulated H/3 width values from the SDSS 
database. Comparison with the final measured values shows that they are not a reliable 
indicator of what is conventionally used as the FWHM of this line. Only about 50% of the 
database values are within 20% of the measured widths. The discrepancies arise from (1) the 
inclusion of the narrow component in the automatic measurement, (2) the departure of the 
profile from a Gaussian, and (3) effects due to contamination by other lines - particularly Fe 
II. 

Aside from the [O III] velocity values that characterize each sample, the statistical 
difference among the samples are limited to a few properties. Among the B, N, and R 
samples, there are no significant differences in the distribution of properties measured. The 
VB sample is differentiated from the other three samples in having narrower H/3 lines and 
broader [O III] A5007 lines. The narrower H/3 lines result in the VB sample having smaller 
black hole mass and larger L/L^d on average than any of the other samples. It is worth 
noting that the VB sample does not differ from the other samples in the distribution of [O 
III] line asymmetry, A20. 
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4. Discussion 

It was previously known that among the objects with blueshifted [O III], Narrow Line 
Seyfert Is (NLSls), which have H/3 FWHM values less than 2000 km s" 1 and are thought to 
be high L/L Edd objects, are overrepresented. Zamanov et al. (2002) studied a small sample 
of objects with extreme [O III] blueshifts, greater than 250 km s -1 relative to the peak of 
the H/3 line, drawn from a heterogeneous sample of 216 objects. These authors found that 
the seven extreme objects all had H/3 FWHM values below 4000 km s" 1 , with three of the 
extreme objects falling in the NLS1 regime. About half of the sample from which these 
objects were drawn has H/3 FWHM values below 4000 km s _1 . 

Both the derived black hole mass and the Eddington ratio (L/L E(M ) depend on the H/3 
FWHM and the luminosity. Limiting the comparison to just the VB and N subsamples, the 
average luminosities of the two are not significantly different but the H/3 widths are quite 
different - 2175 km s _1 for the VB sample vs 4456 km s _1 for the N sample. This leads to 
black hole masses substantially smaller and Eddington ratios substantially larger on average 
for the VB objects. A comparison of the distributions, however, shows that the difference in 
Eddington ratio is far more significant than that in black hole mass (0.8% vs 4.5% probability 
from a student t-test). 

Table 1 also shows, however, that not all narrow H/3, high Eddington ratio objects have 
large [O III] blueshifts and that not all the objects with large [O III] blueshifts have narrow 
H/3 and high Eddington ratios. In order to investigate whether other properties are related 
to the distinction between VB and non-VB objects, another subsample, consisting of 12 
objects that appeared to be NLSls with [O III] velocities near zero, was defined. The same 
H/3 and [O III] characteristics were measured for this sample, called the N-NLS1 sample, and 
these values are also given in Table 1. Other than [O III] velocity, this N-NLS1 subsample 
has properties indistinguishable from the VB sample, except for the width of the [O III] 
line, which has an average of 335 km s _1 , consistent with the N value of 328 km s _1 , and 
inconsistent with the VB value of 523 km s _1 . This is in contrast to the smaller average 
black hole masses calculated for the VB (and N-NLS1) objects, and therefore in conflict with 
the idea that the [O III] width is indicative of the stellar velocity dispersion, which is, in 
turn, correlated with the black hole mass. Thus, one could characterize the situation with 
two statements: (1) objects with higher Eddington ratio are more likely to have large [O 
III] blueshifts; and (2) objects with large [O III] blueshifts have anomalously broad [O III] 
emission lines. 

This result is shown graphically in Figure 4, in which black hole mass is plotted against 
the log of the [O III] line width (a) for the separate subsamples. Each subsample is plotted 
as a different symbol, and the mean values for each subsample are plotted as large symbols 
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of the same type. The line in Figure 4 is the Tremaine et al. (2002) relation between black 
hole mass and stellar velocity dispersion. It is apparent that the VB subsample (solid circles) 
have similar black hole masses to the N-NLS1 subsample (open circles), but that the [O III] 
widths are much larger for the VB objects, which lie substantially below or to the right of 
the line. It is interesting that the N-NLS1 subsample also lies slightly below or to the right 
of the line, suggesting the possibility of a similar, though smaller discrepancy between stellar 
velocity dispersion and [O III] line width. 

It has been claimed that objects with large [O III] blueshifts also are more likely to show 
blue asymmetries in that line. That does not appear to be the case in this sample. There 
is no statistically significant difference in the distribution of A20 values between any of the 
subsamples. Only six of the 59 A20 values do not show a blue asymmetry, and, of these, 
two indicate no asymmetry in either direction. It appears that blue asymmetries are equally 
dominant in objects of all types, both with and without [O III] blueshifts. 

The distribution of equivalent width of the [O III] A5007 line among subsamples is 
also worth noting. While there is no formal difference between the distributions within, 
admittedly small, subsamples, five of the 12 objects in the N sample have equivalent widths 
larger than the largest object in the VB sample. Thus, the N sample is characterized by 
a much broader distribution (standard deviation = 31 A) of A5007 equivalent widths than 
is the VB sample (standard deviation = 6 A). The distributions differ formally at the 90% 
significance level. 

Many previous studies have speculated on the cause of [O III] emission with an extended 
blue wing, a broad underlying blueshifted component, or a blueshifted line center. The 
explanations include both outflow models and inflow models. Outflow models are motivated 
by the idea that objects such as NLSls have characteristics that would support radiation 
pressure as a mechanism for accelerating such flows. Imaging observations (Ruiz et al. 2005), 
which show the [O III] emission coming from a biconical region in many cases, support this 
idea. These models must posit some obscuring screen, depending on the assumed scale of 
the outflow, that blocks the observers view of the material on the far side. This could be 
the 'dusty torus' or it could be a larger, central plane of interstellar material in the host 
galaxy. Inflow models have the advantage that they produce the line asymmetry by mixing 
the obscuring material with the NLR clouds, thus, no additional structure that happens to 
match the size of the outflowing, emitting material is needed. It is the material on the near 
side that is invisible, since the backsides of these clouds are being illuminated. These models 
are supported by the idea that the emitting clouds in NLSls are the most dense - as deduced 
from studies of lines arising from ions such as Fe II, Ca II, and Na I (Ferland and Persson 
1989). Thus, both of these types of models are consistent with the increasing frequency of 
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[O III] blueshifts with increasing L/L^d- 

Although the new findings (as well as previous investigations) do not clearly distinguish 
between the two models, a number of interesting conclusions may be drawn. Radial motion 
of narrow line emitting clouds is present in almost all objects. This is the only way to account 
for the common occurrence of line peak blueshifts and the ubiquitous presence of the line 
asymmetry. In the inflow model, we never see the near side of the clouds on the near side 
of the nucleus; in the outflow model, we never see the clouds on the far side of the nucleus. 
This latter point is interesting because it implies that the scale of narrow line emission is 
not large relative to some obscuring plane in the center of the active galaxy, whether that 
is a structure related to the accretion flow or it is the central plane of the galaxy. This also 
suggests that the [O III] width is not generally a suitable surrogate for the stellar velocity 
dispersion, but rather, represents the radial velocity gradient in the outflow or inflow. 

Why have previous studies found the [O III] line width to be correlated with stellar 
velocity dispersion? One possible explanation is that the outflow (or inflow) of the narrow- 
line emitting material is affected by the black hole mass, either through gravity or through 
radiation pressure. If the [O III] line width is indicative of the maximum radial velocity, 
then it would be correlated with the black hole mass. At the same time, the stellar velocity 
dispersion is correlated with the black hole mass, and the two widths appear to be correlated 
with each other. 

The most apparent flaw in this picture comes from the objects with extreme [O III] 
blueshifts. These are predominantly objects with small black hole masses and yet they have 
very broad [O III] lines. Furthermore, they are, in all other ways, identical to another group 
of objects, the N-NLS1 subsample. Clearly another parameter is required. One possibility 
for this parameter is orientation; objects with larger blueshifts and broader lines are those 
in which the outflow is directed closer to the line of sight. 

Putting these pieces together yields the following picture: At a given L/L^dd, outflow 
(or inflow) velocity is larger in objects with larger black hole mass. At a given black hole 
mass, outflow velocity is larger in objects with larger L/Lg^. In addition, as one goes to 
larger L/L Edd , the cone angle of the outflow decreases, and in the objects with large L/L Edd , 
one sees a range from small outflow velocities (in objects in which the axis of the outflow is in 
the plane of the sky) to large outflow velocities (in which the axis is along the line of sight). 
Whether this picture extends all the way to truly edge-on objects in which one sees both 
sides of the outflow, and consequently symmetric lines with twice the flux, is not clear given 
the small number of objects studied here. This model implies that the other parameters that 
distinguish objects - H/3 width, Fe II strength, [O III] equivalent width - are not primarily 
orientation dependent. 
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An interesting comparison can be made with the velocity distribution of the peak of 
the C IV line, studied in higher redshift SDSS quasars by Richards et al. (2002). Working 
with composite spectra generated from large samples of objects divided by C IV blueshift 
relative to the Mg II line, these authors found that more blueshifted C IV lines occurred as 
a result of the red wing of the line being obscured. They suggested that this obscuration is 
due to a central screen of finite geometrical thickness, which is optically thin seen face-on 
but becomes thicker with increasing angle. Note that this is, in some sense, opposite to the 
scheme proposed here, in which the objects with the largest blueshifts are those seen closest 
to face-on. In any case, the geometry and kinematics of the C IV-emitting clouds and the 
[O III] -emitting clouds may be different, and so the explanations may not be inconsistent. 



5. Conclusions 

Characteristics of AGNs with blueshifted [O III] lines have been reported previously. 
However, conclusions have been based on small samples, with heterogeneous data, or requir- 
ing dubious assumptions. In this study, a large, uniform sample of AGNs from the SDSS 
DR1 have been analyzed. Objects with discrepant [O III] redshifts have been identified 
and subsamples constructed to explore the properties of objects populating different regions 
in the [O III] relative velocity distribution. Analysis of these samples yields the following 
conclusions. 

1) A significant fraction, as many as half, of all AGNs show [O III] emission that is 
blueshifted to the point that the line peak is shifted by tens to hundreds of km s _1 . 

2) This shift is a true shift relative to the systemic velocity, as determined from the low 
ionization narrow forbidden lines. 

3) The AGNs with the largest blueshifted [O III] lines tend to have high L/L Edd: using 
black hole masses determined from their luminosities and H./3 widths. However, this is not 
a tight correlation; there are high L/L Edd objects that do not show blueshifts, and lower 
L/Ls^ objects that do. 

4) The single anomalous property common to all the objects with highly blueshifted [O 
III] emission is the large width of the [O III] lines. These average more than 50% broader 
than those in similar objects that have [O III] close to the systemic velocity. 

5) Conversely, the blueshifted [O III] lines do not tend to be any more asymmetric or 
any weaker in strength than the non-blueshifted ones. All the subsamples show a similar 
tendency to have a sharp red edge to the line and a blue edge that is shallow, or has a wing 
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or a shoulder. 

While these characteristics do not identify a unique model, simple arguments suggest 
that the maximum radial velocity of the [O III] emitting clouds is controlled by both black 
hole mass and Eddington ratio. A radial velocity gradient is primarily responsible for the line 
width, which then allows an orientation parameter to play a role. The apparent correlation 
of [O III] line width with stellar velocity dispersion in previous studies may be due to the 
fact that both of these widths are somewhat dependent on black hole mass. 
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Fig. 1. — Distribution of absolute magnitude in the I band (top) and redshift (bottom) for 
the sample of 804 SDSS DR1 quasars studied. The dashed line indicates those for which any 
low ionization lines were measured. The dotted line indicates those for which lines from two 
or more low ionization species were measured. 
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Fig. 2. — For each of the three ions indicated, the distribution of redshift relative to that of 
other low ionization lines. 
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Fig. 3. — Relative redshift distributions for all measurements of low ionization lines (bottom) 
and for the [O III] A5007 line (top). The dotted line in the top panel shows the objects for 
which two or more low ionization species are used to calculate the systemic velocity. 
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Fig. 4. — Black hole mass plotted against the [0 III] line width. The VB subsample is plotted 
as solid circles, the B subsample is plotted as open squares, the N subsample is plotted as 
solid triangles, the R subsample is plotted as crosses, and the N-NLS1 subsample is plotted 
as open circles. The large symbols of each type show the position of the mean values for 
that subsample. The line is the Tremaine et al. (2002) fit to the relation between black 
hole mass and stellar velocity dispersion. 
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6 


38 
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19 
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6 


41 
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